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BERNARD SHIFFMAN AND STEVE ZELDITCH 

Abstract. We introduce several notions of 'random fewnomials', i.e. random polynomials 
with a fixed number / of monomials of degree N. The / exponents are chosen at random 
and then the coefficients are chosen to be Gaussian random, mainly from the SU(m + 1) 
ensemble. The results give limiting formulas as N — > oo for the expected distribution of 
complex zeros of a system of k random fewnomials in m variables (k < m). When k = m, for 
SU(m+l) polynomials, the limit is the Monge- Ampere measure of a toric Kahler potential on 
CP™ obtained by averaging a 'discrete Legendre transform' of the Fubini-Study symplectic 
potential at / points of the unit simplex £ C R m . 
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Introduction 

This article is concerned with the distribution of complex zeros of random systems of 
fewnomials. Fewnomials are polynomials 

£ c a z a : #{a:c a ^0} = f «N, 

a:\a\<N 

with a relatively small number f « N of non-zero coefficients in comparison to the degree. 
For instance, z 10 > 000 + 5Z 1 ' 000 - 3z 100 - 1 is a fewnomial but z 10 > 000 + z 9 < 999 + ■■■ + z + 1 
is not. The fundamental idea is that the number of monomials, rather than the degree, 
measures the complexity of a polynomial system [Kh]. The purpose of this article to begin 
an investigation of fewnomial complexity bounds from a probabilistic viewpoint. In this 
article, we introduce several natural ensembles of 'random fewnomial systems', and study 
the expected distribution of their complex zeros. In subsequent articles, we plan to study 
real zeros of real fewnomials and the more difficult problem of the correlations and variance 
of both real and complex zeros. The overall purpose is to study from a statistical point 
of view Khovanskii's bounds [Kh] on the Betti numbers of real algebraic varieties given by 
fewnomials and on the number of real zeros of a full fewnomial system. Statistical properties 
of the topology of real algebraic varieties given by non-fewnomial polynomials (i.e., with all 
the coefficients nonzero) have been studied, for example, in [Bu, Ro, ShSm, GW] and the 
references in these papers, but not much is known statistically about zeros of fewnomials. 

To put our problem into context, let us recall Khovanskii's theorem: Let P = (Pi, ... , P m ) 
denote a system of m complex polynomials on (C*) m , and let Aj = Ap j denote the Newton 
polytope of Pj, i.e the convex hull of the exponents appearing non-trivially in Pj. Let U C T m 
be an open set, where T m C (C*) m is the real m-torus, and let N(P, U) be the number of 
zeros with arguments lying in U. When U = T m , N(P, U) counts the total number of zeros 
in (C*) m , which by the Bernstein-Kouchnirenko theorem [Be, Ko] can be expressed in terms 
of the mixed volume V(Ai, . . . , A m ). Given an angular sector determined by U, the number 

S(P, U) := V(A 1 , A m )Vol(f/)/Vol(T m ) (1) 

may be viewed as the 'average number' of complex zeros in the sector among random polyno- 
mial systems with the prescribed Newton polytopes Aj. We denote this class of polynomial 

systems by 

P e Poly(A 1 , . . . , A m ) := {(P u . . . , P rn ) : A Pj = A,}. 
Khovanskii's complex fewnomials theorem [Kh, §3.13, Th. 2] asserts that 

sup \N(P, U) - S(P, U) | < U(U, Ax, ... , A m )<p(m, /) (2) 

PePoly(Ai,...,A 

where / is the number of non-zero coefficients of the system, and where H(U, Ai, . . . , A m ) 
is the smallest number of translates of a certain region A* C T m required to cover the 
boundary of U . One of the principal applications of this result is to give an upper bound 
for the number N^(P) of real zeros of a fewnomial system: If Uj is a sequence of small balls 
around {0} shrinking to the point {0}, one has II = 1 and S(P, Uj) — > and one obtains a 
bound of the form 

\N R (P)\ < V (mJ) (3) 
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entirely in terms of the number of non-zero monomials appearing in it and not its degree. 
We will refer to / as the fewnomial number of the system. 

Khovanskii's result may be interpreted in terms of the angular projection Arg : (zi, . . . , z m ) - 
(j^j, . . . , -^j) of the zero set to the real torus T m . His result (in the full system case) says 
that the angular projection of the fewnomial zero set is rather evenly distributed in T m . As 
a result, not too many zeros concentrate on the real set where 6 = 0. Note that his measure 
of the concentration, taking the supremum in (2), is very astringent and is governed by the 
extreme cases. The idea of our work is to study its average value over fewnomial systems 
and polynomials. 

The motivation for the statistical study is that the known estimates of (p(m, /) are very 
large and are widely conjectured to overestimate the the bound by many orders of magnitude. 
Khovanskii's bound states that <p(m, f) < 2 m 2^~ 1 )/ 2 (m + 1)'. See [BBS, BRS] for relatively 
recent bounds and [Sot, Stl, St2] for further background. A conjecture of Kouchnirenko, 
as corrected and refined by a number of people, states that the maximum number of real 
zeros in the positive real quadrant should be roughly \f\ 2m where |/| is the total number 
of monomials in the system. The uncertainty as to the true order of magnitude of ip(m, /) 
suggests studying the bound probabilistically. The bound (3) resembles a variance estimate 
although it is measured in the much more difficult sup norm. It reflects the extremal behavior, 
which may only occur very rarely. This raises the question, what is the expected or average 
order of magnitude of the variance? 

In this article we begin the study of random fewnomial systems by introducing several 
probability measures on spaces of complex fewnomials — i.e., on the set of pairs (S,P) of 
spectra and polynomial systems with the given spectra. Our main results give the expected 
limit distribution of complex zeros in the ensembles. For example, Theorem 4 says that 
for a random system (Pf, . . . , P^) of fewnomials on C m , each of fewnomial number / and 
of degree N, where the exponents are chosen uniformly at random and the coefficients are 
chosen at random from the SU(m + l) ensemble (described below), the expected distribution 
of zeros in (C*) m is asymptotic to 

(a2 r r -, \ dO df) 

—— / max (p, A'"> - (AUog A' > dX 1 ■ ■ ■ dXH dp, ■ ■ ■ dp m -± ■ ■ ■ - 

where z = (e pi / 2+%ei , . . . ^ e Pm/2+id m ^ £ j s the unit simplex in W 71 with probability measure 
dX = ml dXi ■ ■ ■ dX m , and A = (1 — |A|, Ai, . . . , A m ). 

0.1. Fewnomial ensembles. We consider several natural definitions which are motivated 
by different kinds of applications. More precise and detailed definitions are given in §2. 

We denote the space of all complex holomorphic polynomials of degree N by Poly (AT). 
By the spectrum (or support) of a polynomial P, we mean the set Sp of exponents of its 
non-zero monomials. We denote the space of polynomials with spectrum contained in S by 

Poly(S) = {P(z u ...,z m ) = Y,Cc l Xa(z), X a(z):=z a }, ScN m . (4) 

The Newton polytope of P is the convex hull A P of the spectrum Sp. More generally, we 
consider a system of k < m polynomials Pi, . . . , P^ mm complex variables, and write 

Po\y(S 1 ,...,S k ) = {(P 1 ,...,P k ): P 3 e Poly(^)}. (5) 
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When k = m we speak of a 'full' system, where the simultaneous zeros are almost always 
zero-dimensional . 

In all of our definitions of random fewnomial system, the numbers fj of elements of each 
spectrum Sj and the degrees N of the polynomials Pj are fixed. We then randomize with 
respect to the spectra Sj and with respect to the coefficients Cj a . With regard to the spectra, 
there are several natural choices of probability measure: 

(I) Fixed spectrum up to dilation: Here, we fix a spectrum S, and then dilate it determin- 
istically as the degree N grows, i.e. scale S — > NS. This notion of random fewnomial 
is analogous to our notion of random polynomial with fixed Newton polytope in [SZ2]. 
The main difference is that the 'polytope' which we dilate is non-convex; indeed, it 
just consists of a fixed set of / points. In fact, the techniques of [SZ2] generalize quite 
naturally to all non-convex polytopes. The only randomness is then with respect to 
the coefficients. The result is given in Theorem 1. 
(II) Dilates of a random spectrum from a polytope A: in this ensemble, we fix A, choose 
the spectrum at random from A and then dilate the resulting spectrum. It is only 
a small step from case (I), but is apparently important in computational work. (We 
thank Maurice Rojas for emphasizing the interest of this case.) The result is given 
in Corollary 3. 

(III) Random spectra of degree N: At the opposite extreme, we may choose the spectra 
completely randomly (with respect to counting measure) from all possible /-element 
sets of exponents a G N m of length \a\ = \ai\ + • • • + \a m \ < N — i.e., subsets of 
the integral simplex Z m H NT,, where T = {x G M. m : Xj > 0,Yl x j — 1} * s the 
unit m-simplex, and iVE is its dilate by N. We put uniform measure on Z m D NT 
and then choose spectra Sj C Z m fl NT of fixed cardinality / independently and 
uniformly relative to counting measure. The result is given in Theorem 4 for SU(m + 
1) fewnomials and in Theorem 7 for general toric Kahler potentials. 

(IV) Random spectra contained in fixed Newton polytopes: Rather than just consider 
the simplex, and motivated by Khovanskii's variation formula, we now fix k convex 
lattice polytopes A 1; . . . , A k and then choose random spectra Sj C Aj fl Z m with 
fixed cardinalities fj independently with uniform measures from these polytopes. We 
then replace the Aj by their dilations by NAj to obtain higher degree ensembles. 
We note that the convex hull of Sj is contained in Aj but equality rarely occurs. The 
result is given in Theorem 6. 

Having decided on an ensemble of spectra a, we then define probability measures on the 
coefficients c a . We only consider Gaussian probability measures and make standard choices 
which are consistent with Khovanskii's bound. A key point is that Gaussian measures are 
determined by inner products in the space of polynomials. We choose the inner products as 
in [SZ1, SZ3, SoZ2, SoZ3] to be those G N (<p,v) of the weighted L 2 spaces L 2 (CF m , e~ Nv du) 
of pluri-potential theory, which are specified by a Kahler potential ip or Hermitian metric 
h = e _</3 and a measure dv on CP m . It is natural to restrict to tp, v which are toric, i.e. 
invariant under the standard T m torus action on CP" 1 . Then the monomials {z a } are always 
orthogonal and the Gaussian ensembles only differ in the L 2 - norms 

Q GN{ ^) = m\l N ^)= [ \z a \ 2 e- N ^du(z) (6) 
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of the monomials, viewed as homogeneous polynomials of degree N (so that \a\ < N); 
equivalently, the Gaussian measures only differ in the variances of the coefficients in the 
monomial basis. We refer to §1 for details. 

Since our emphasis is on the fewnomial aspects we only consider some basic examples of 
((p,v). In particular, we concentrate on the model case of SU(m + 1) polynomials, where 
(p(z) = log(l + ||-2|| 2 ) is the Fubini-Study potential and where dv = ^j(^dd(p) m is the 
Fubini-Study volume form. 

Given the inner product Gn((p,v) underlying the Gaussian measure, we normalize the 
monomials to have L 2 -norm equal to one, by putting 



QG N {<p,v){Oi) 

and then express polynomials of degree N as the orthonormal sums 

The Gaussian measure 7 at induced by Gn(<£, v) is defined by the condition that the c a are 
independent complex normal variables of mean zero and variance one. 

The Gaussian measure djN on Poly(iV) induces conditional Gaussian measures 7^5 on 
the spaces Poly (5); i.e., 

dj N \s(PN) = ^ e ~ E|cQ ' 2 dc, P N = J2 c ^a, (9) 

aeS 

where f a is given by (7). Probabilities relative to ^n\s can be considered as conditional 
probabilities; i.e., for any event E, 

Prob 7jv {P eE\S P = S} = Prob 7jv|S (£). 

We denote by Hjn\s the expectation with respect to the conditional Gaussian measure 7jv|s- 
For a further discussion of conditional probabilities on polynomial (and more general) en- 
sembles, see [SZZ] 

Some of the possible (and well-studied) choices of the inner products and Gaussian mea- 
sures are the following: 

(a) The SU(m + 1) ensembles defined above. On all of Poly(iV) the expected distribution 
of zeros for each N is uniform with respect to the SU(m + l)-invariant volume form 
on CP m (i.e. the Fubini-Study form). 

(b) General toric Gaussian measures induced by T m -invariant Hermitian metrics h = 

on the line bundle 0(1) — > CP m with positive curvature form = iddtp (i.e. with a 
plurisubharmonic weight) and with v = dV v := ■^j(^dd(p) m . We suppress geometric 
notions in this article, but state the general result in Theorem 7 (see §4). 

(c) The m-dimensional Kac-Hammersley ensembles, where Qn{®) = 1- Here, tp = 
(the opposite extreme from subharmonic weights) and v = 5x m - The norms of the 
monomials are independent of N and only involve a fixed inner product on C m . In 
dimension one, the zeros of degree N polynomials (with full spectrum Z n [0,7V]) 
concentrate on the unit circle as N — > 00 [Ha], and in dimension m the zeros of 
degree N polynomials concentrate on the torus T m [BS]. We briefly discuss this 
ensemble in 85. 
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Remark: Khovanskii [Kh] and Kazarnovskii [Kal, Ka2] consider ensembles where one fixes 
the spectra Sj and chooses coefficients at random from the ensemble 

CplSil- 1 x CP |52hl x • • • x CP |5m|- \ (10) 

i.e., the product projective space of coefficients of polynomials with the prescribed spec- 
tra, equipped with the probability measure obtained by taking the product of (normalized) 
Fubini-Study volume measures on the factors. These ensembles amount to choosing the 
complex coefficients at random from the Euclidean spheres S* 2 ' 53 ' -1 and are easily seen to be 
equivalent to Gaussian random polynomials J2 a CaZa with c a independent complex normal 
variables of mean zero and variance one — i.e., they are equivalent to the Kac-Hammersley 
ensembles described above. 



0.2. Expected distribution of zeros. Having fixed an ensemble of fewnomials, our inter- 
est is in the configuration of zeros 

Z Pl ,..., Pk := {z e (CT : Pi(z) = ■■■ = P k (z) = 0} 

of a random fewnomial system with k < m. Here, C* = C \ {0}. We refer to k = 1 as the 
random fewnomial hypersurface case and to k = m as the point case. 

To each zero set we associate the current of integration [Z Pli ,„ jPk ] e V' k ' k ((C*) m ) over the 
zeros of the system: 

([z Pl _ Pk ] , = / l> , ^ g v m - k > m - k ((c*) m ) . 

Jz Pl ,...,p k 

In the point case, [Z Plj ... 5 p m ] is obtained by putting point masses at each zero, 

[Z Pl ,..., P J = S * ' 

zeZp!,...,p m 

and the expected distribution is determined by the expected values of the random variables 

M%(P?, -,PZ)--= [Z Pi ,.,pJ (U) = #{^f/: P»(z) = ■ ■ ■ = P%(z) = 0} 

counting the number of zeros in an open set U C (C*) m . 

The expected distribution of zeros varies widely among the ensembles above. This is not 
surprising if one recalls, for instance, that zeros of random Kac polynomials concentrate on 
the unit circle, while those of SU(2) polyonomials are uniform with respect to the standard 
area form of CP 1 = C U oo, while those of polynomials with fixed Newton polytope have 
a forbidden region where zeros have an exotic concentration. In particular, the 'average 
number' S(P, U) of zeros in the angular sector U considered in Khovanskii's variance estimate 
(2) is itself a random variable which depends on the convex hull of the spectrum of P. 

0.3. Statement of results. We will consider the zero distribution as a measure on (C*) m = 
IR™ x T m and denote points by z = e p l 2+%e in multi-index notation. Here, T m denotes the 
real torus T m = {S 1 ) rn C (C*) m . Given a locally bounded plurisubharmonic function (p we 
denote by MA {(p) the associated Monge- Ampere measure 



(\ m 
^ddp) G p' m - m (CP m ). 
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When ip is invariant under the T m action on (C*) m , then 



MA(<p) = det^—D 2 p <pjdpde, (11) 

where D 2 is the real Hessian on M. m . 

Our results are asymptotic formulas as the degree iV — > oo, but with the number / of 
monomials held fixed. For each ensemble, the limit distribution of zeros in the point case 
is the Monge-Ampere measure of a limit T m -independent potential, and thus the formula 
is of the type (11). The results are very similar for the Fubini-Study SU(m + 1) ensemble 
and for general toric Gaussian measures based on inner products Gn(<p, dV^) with p a toric 
Kahler potential. Hence we concentrate on the SU(m + 1) case, and only briefly indicate the 
modifications needed for the general toric Kahler case. 

Our first result concerns the ensemble with dilates of a fixed spectrum. Since the lattice 
points lie in M. m we use upper subscripts to index the different points in the spectrum and 
lower subscripts to index their coordinates. We recall that Eat|s refers to the expectation 
with respect to the conditional Gaussian measure Jn\s- 

Theorem 1. Let S = {A 1 , . . . , A-^} be a fixed spectrum consisting of f lattice points in pE. 
For random m-tuples (Pf,...,P^) of fewnomials in Po\j(NS) , with coefficients chosen 
from the SU(m + 1) ensembles of degree pN , the expected distribution of zeros in (C*) m has 
the asymptotics 

N- m E NplNS [Z PiN _ Pj ,] ^ p m MA (max [(p, A) - (A p ,log>)]) . 

Here, X p = (p - |A|, A i? . . . , A m ) and log X p = (log(p - |A|), log A i? . . . , log A m ). 

For a spectrum S C pE, we let C p s denote the Monge-Ampere potential in Theorem 1: 

£* 8 (p) := max Up, A) - (>, log A>>1 , p G M m . (12) 

It is kind of discrete Legendre transform of the entropy function (A p , logA p ), which is the 
symplectic potential corresponding to the Fubini-Study Kahler potential. 

We note that the expected limit distribution is a singular measure invariant under rotations 
of the angular variables and supported along the O-dimensional corner set of the piecewise 
linear function C p s (p). This reflects the heuristic principle that the zeros of a fewnomial 
should come from its sub-fewnomials with fewnomial number / = m + 1 . 

With no additional effort, we could fix the spectra separately for each polynomial in the 
system, and obtain: 

Theorem 2. Let S 1 , . . . , S k be fixed finite spectra consisting of lattice points in pE 7 where 
1 < k < m. For random fewnomial k-tuples (Pf, . . . , P£) in Poly (NS 1 ) x • • • x Poly(A^S' fc ) ; 
with coefficients chosen from the SU(m + l) ensembles of degree pN , the expected zero current 
in (C*) m has the asymptotics 

N- k E NSl _ NSk [Z PiN _ PkN ] /\ (^ddC%{p)\ . 

T = 1 ^ ' 
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We now state the result for fewnomial ensembles in which we randomize the spectra in 
the sense of (II): 

Corollary 3. Let A C j)S be a (fixed) Newton polytope, let S 1 , . . . , S k be random spectra 
contained in A with fewnomial number f , and let P^ , . . . , P^ be random fewnomial k-tuples 
(Pf, P£) in Poly(iVS' 1 ) x • • • x Po\j(NS k ), with coefficients chosen from the SU(m + l) 
ensembles of degree pN. Then the expected zero current in (C*) m has the asymptotics 



y ' J ' \ seC(Aj) 



Next, instead of dilating random spectra, we consider completely random spectra as de- 
scribed in (III) and we obtain: 

Theorem 4. Let 1 < k < m, and let (Pi, . . . , P k ) be a random system of fewnomials of 
fewnomial number f and of degree N, where the spectra Sj are chosen uniformly at random 
from the simplex NT, and the coefficients are chosen from the SU(m + 1) ensemble. Then 
the expected zero current in (C*) m has the asymptotics 

N- k E Nf [Z pN P N}^(^dd[ max \(p, \ j ) - (V, log V)] d\ l ■ ■ ■ d\ f \ . 

1 k \2n Jxf j=i,.. .J L J J 

Here, A = A 1 = (1 — |A|, Ai, . . . , \ m ), d\ = m\ d\\ • • • d\ m . 

The limit measure is thus the Monge- Ampere measure of the limit potential obtained by 
averaging the discrete Legendre transform £| Ai \ f y(p) from Theorems 1 and 2 (with p — 1) 

over all choices of points A 1 , . . . , A^ of E. 

We note that the averaging smooths out the corners. Indeed, we have the following more 
explicit formula for the expected limit distribution: 

Corollary 5. Let (P 1 , . . . , P k ) be the random system of Theorem 4- Then 

N- k E NJ [Z PiN _ PkN ] | WFS -±-J~dd([l- D b (t; p)] f )dt} k , 



where 



D b (t; p) = ml Vol ({AGS: (A, log A) - (p, A) + log (1 + \e p \) < t} 



Here, u FS = ^dd\og(l + \e p \) is the Fubini-Study Kahler form on (C*) m C CP m . The 
quantity D b {t] p) is the distribution function for the pointwise logarithmic decay rate b\(p) 
of the monomials ifNx (see §2.3), regarded as a random variable (with parameter p) on E. 
Note that the integral in Corollary 5 is actually over a bounded interval. 

We can also generalize Theorem 4 to the ensemble (IV): 

Theorem 6. Let 1 < k < m, let Ai C piE, . . . , A k C p k ^ be Newton polytopes and let 
(P±, . . . , P k ) be a random system of fewnomials of fewnomial numbers fi, ■ ■ ■ , fk respectively, 
where the spectra Sj are chosen uniformly at random from the simplices NAj and the coef- 
ficients are chosen from the SU(m + 1) ensemble. Then the expected zero current in (C*) m 
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dX l ---dX f * 



N- m E[Z pN pN }^ A ( p- -. dd [ max \{p,X l ) - (X 1 ,log A' 

L Pl >">*>> 1 /J y 2,71 V0l(Aj)^ '=!.-./* I- 

i/ere, A = A 1 = (1 — |A|, Ai, . . . , X rn ), dX = dXi ■ ■ ■ dX m . 

The key analytical ingredient in the proofs of these results is an asymptotic formula for 
the expected mass density of the above systems of random polynomials as N — > oo. It 
is given by the conditional Szego kernel IIjv,q|s with respect to the norm Qa N (<p,v) (6) and 
spectrum S, i.e. the kernel of the orthogonal projection (Szego kernel) onto the subspace of 
polynomials under consideration: 

Vn, qis (\p(z)\ n ,) = Eni^ = n w (^) . (13) 

Thus, the results depend on the asymptotics of the Szego kernels Un,q\s(z, z). 



0.4. More general toric weights. We briefly indicate the generalization when the SU(m + 
1) (Fubini-Study) inner product on Poly(A^) is replaced by Gn^, dV^) for a general toric 
Kahler potential ip. 

The polytope P of the toric variety is defined by a set of linear inequalities 

£ r (x) := (x, v r ) — A r > 0, r — 1, d, 

where v r is a primitive element of the lattice and inward-pointing normal to the r-th (m — 1)- 
dimensional facet F r = {£ r = 0} of P. 

A T m -invariant Kahler potential on (C*) m defines a real convex function on p E M. m . Its 
Legendre transform 

u v (x) := C<p(x) := sup ((x, p) - p(e p )) 
p 

is the symplectic potential u v . Equivalently for x G P, there is a unique p such that 
Vpip = x, and u v (x) = (x,p x ) — (p(p x )- In the Fubini-Study case, P — E, ip — log(l + e p ), 
and ufs(x) — Sfc^fc( x ) l°g^fc( a; ) where £k(x) — Xk for k — 1, . . . , m and £ m +i(x) — 1 — \x\ 
where \x\ — x\ + • • • + x m (in multi-index notation on W 11 ). Thus, 

ufs(X) = (A, log A). 

The Kahler potential is the Legendre transform Cu v (p) of its symplectic potential. If we 
allowed all possible spectra in the ensemble (hence not a fewnomial ensemble), the discrete 
Legendre transforms with respect to /-element subsets would converge to the usual Legendre 
transform and the potential in (12) would become ip. Thus, the impact of the restriction 
to / monomials is that in place of the Legendre transform we have an average of discrete 
Legendre transforms. 

As this indicates, the result for a general toric Kahler Gaussian ensemble for CP m and 
polytope E, defined by Gj^(p, dV^), is the following: 
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Theorem 7. Consider the ensembles of type (III) as in Theorem 4, but with Gaussian 
measures induced by the inner product Gn{<£, dV v ) corresponding to a toric Kahler potential 
on CP m . Then the expected distribution of zeros in (C*) m has the asymptotics 

N- m E[Z pN pN ] ->• l—dd [ max Up, X J ) - uJX)] dX 1 ■ ■ ■ d\ f 

1 ' ' k X^TT Jy: I I / 

The proof is almost the same as for the Fubini-Study case and is indicated in §4. In §5, 
we also indicate the modifications in the case of the fewnomial Kac-Hammersley ensemble. 

1. Preliminaries 

In this section, we review the relation between inner products on spaces of polynomials 
and associated Gaussian measures on the space. The inner products implicitly involve a 
choice of Kahler metric on CP m . The associated Kahler potential determines the shape of 
the modulus of each monomial and its concentration properties. 

We may identify a not necessarily homogeneous polynomial / on C m of degree < N by its 
homogenization polynomial 

F((o, •..,&»)=£ CaC A (C A = Co° • • • ct) 

\\\=N 

of degree N in m + 1 variables, where 

f(z 1 ,...,z m )=F(l,z 1 ,...,z m ) = c a z a (z a = z?---z%»), 

\a\<N 

where c a = C & n, a N = (N — \a\, a±, . . . , a m ), \a\ — Y^j=i a j- Homogeneous polynomials of 
degree iV on C m+1 are equivalent to holomorphic sections H°(CF m , O(N)) of the Nth power 
of the hyperplane section bundle. This geometric identification is useful in interpreting the 
concentration properties of monomials in terms of curvature. 

We let e G H°(CF m O(l)) be the degree 1 polynomial e (Co, • • • , Cm) = Co- Then e is a 
local frame over the affme chart Uq = {Co 7^ 0} ~ C m . We fix a Hermitian metric h on 0(1). 
In the local frame eo, the metric has the local expression h = e~ v , where tp is known as the 
Kahler potential. The Kahler form is denoted by u v = ^ ddip. 

We define the inner product on Poly(iVS): 

(/, 9)h = —yf f{zjgJz)e- N ^ cu™(z), f, g G Poly(TVE). (14) 

ml J C m v 

The inner product is determined by the matrix of inner products on the distinguished basis 
of monomials Xa- All of our inner products are T m -invariant and hence the monomials are 
automatically orthogonal. The inner products are then determined by the norming constants 
(6), specifically, 

Q(a) = Q GNiv , dVv )(a) = —f \z a \ 2 e- N ^u™(z). (15) 

The inner product induces a Gaussian measure 7^ on any subspace S C Poly(A r S). Again 
assuming that the monomials are orthogonal, the basis (7) is (, )h orthonormal and we may 
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write any polynomial in the form 

P N = ^2 c^Va- 

The associated Gaussian measure is defined by the condition that the coefficients of this 
orthonormal expansion are independent complex normal random variables. 

The Szego kernel (or weighted Bergman kernel) for the line bundle O(N) with metric 
h N = e~ Nlp is given over C m by 

U NjQ (z, w) = e- N " <fi • (16) 

It is the kernel for the orthogonal projection from L 2 (X) — > H (CF m ,O(N)), where X — > 
CP m is the unit circle bundle in (L*,h*) with fibers X z = {e^ +ie (l, z , . . . , z m ) : 6 G R} 
over points z G C m ; see [SZ1]. For spectra S C Z m fl iVE, then the kernel for the orthogonal 
projection L 2 (X) — > Poly(S') is the conditional weighted Bergman kernel given by 

U NtQls (z, w) = e- N " V N a (z)WM • (17) 

1.1. The SU(m+l)-ensembles. This is the Gaussian ensemble defined by the inner product 
arising from the Fubini-Study metric cp = log(l + \z\ 2 ). Then u FS = -£^dd\og(l + \\z\\ 2 ) is 
the Fubini-Study Kahler form on C m C CP m and 



^) = ^( T^ofpW^W, /^GPoly(TVE),. (18) 
mi j£m yi -\- \\z\\ ) 

The norming constants for the inner product (14) are: 



|Xa|| \/(Xco Xa) 



(N + m)\( N a ) 



a J (N — \a\)\ai\ ■ ■ • a r 



(19) 



Thus we have an orthonormal basis for Poly(iVE) given by the monomials 



,1 HN + m)\(N 



In this case, the circle bundle X is the unit sphere S 2m+1 C C m+1 . We now regard 
the sections of H°(CF m , O(N)) as homogeneous polynomials restricted to X = S 2m+1 . By 
identifying the point z G (C*) m with the lift x = n^n^Wu^ (1? z i, ■ ■ ■ > -O £ 5' 2m+1 , we may 
write the homogenized monomials on 5' 2m + 1 i n affine coordinates (zi, . . . , z m ) as 

The corresponding L 2 normalized monomials are then: 



(N + m)\ fN\ z a Ar , , 
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In short, the SU(m + 1) ensemble of random polynomials of degree iV consists of polyno- 
mials of the form 




(N + m)\ (N 



N\ \a (1 + \\z\\ 2 ) N l 2 ' 



(23) 



where c a are independent complex normal variables of mean zero and variance one. 

Specializing (17) to the Fubini-Study metric, we have the following definition (where we 
omit the subscript Q indicating the norming constants): 

Definition: Let S C 7L m fl NT. The conditional Fubini-Study Szego kernel Hn\s is the 
kernel for the orthogonal projection to Poly(S') with respect to the induced Fubini-Study 
inner product: 

U Nls (x, y) = £ j^rX WWW) = E *W • (24) 

The conditional Szego kernel can be written explicitly on C m as 

Un ^ W) = M (1 + ||^||^/2(i + || w ||2)N/2 • ( 25 ) 

It is the two-point function for the conditional Gaussian ensemble Poly (S) C Poly (AT). 
The full Fubini-Study Szego kernel is given by 

(N + m)\ ^ ^ „^ (N + m)\ E| a |<iv&V 

|<*|<JV ' v II II V V II H ) 



(N + m)\ 
M 



1 + (z,w) 



N 



;i + ||^|| 2 ) i /2(i + \\ w \\ 2 y/\ 

2. Fewnomial Ensembles 



(27) 



2.1. Precise definitions of random fewnomials. We now define more precisely the en- 
sembles which allow for any fewnomial system. We fix the degree N, and first consider the 
case of one random fewnomial. We specify a set of lattice points by its characteristic function 

a : NT n N m ->■ {0, 1}, (28) 

which may be regarded as an occupation number, designating whether a lattice point is 
occupied (cr(a) = 1) or unoccupied (cr(a) = 0. We denote by \a\ = XLeJVE a ( a ) the number 
of elements in the set, and by Supper = {a : a(a) = 1} the support of a. We put: 

C NJ = {a : NT n N m -> {0, 1} such that \a\ = /}, (29) 

and we denote the number of such subsets by 

C(NJ) = \C NJ \ =[ {m f >) = j^yyj, N-f + OiN-f- 1 ). (30) 
A polynomial with (at most) / non-zero terms can then be written in the form: 

p*M= W\ = f- (31) 
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Thus the space of random fewnomials is given by: 

T NJ = {(a,P) G C NJ x Poly(iV) : P G Poly(Supp a)}. (32) 

There is a natural projection n : F N j — > C^j taking (a, P) — > a and the 'fiber' of this 
projection is Poly(Supp a). The set of fewnomial systems of m polynomials in m variables 
with fewnomial numbers (/i, . . . , f m ) is then given by 

^N,{h,-,fm) '■= ?N,h x • • • x 3~N,f m ■ (33) 

It is also natural to consider fewnomials with spectra contained in a given Newton poly- 
tope. We therefore fix a convex lattice polytope A C (for some p) and replace £ every- 
where by A. Thus, we define 

Cnja = W : NA n ^ m -> {°> !} such that W\ = /}> ( 34 ) 

and 

•Fjv./.a C C nja x Poly(iV) = {(a, P) : Supp (a) C NA, P E Poly(Supp a)}.. (35) 

Similarly we define F N ^ h ,...j m)A ^.. Am for systems. 

We now induce probability measures on J^nj and Fnja^ °y regarding them as 'fibering' 
over C(N, f), by putting counting measure on Cnj and by putting the conditional measures 
dlN\sup P a on the 'fibers'. 

Definition: The ensemble of random SU(m + 1) fewnomials of degree N and fewnomial 
number / is the space J-'nj endowed with the probability measure dfiNj defined by 

f 9(S,P)d» NJ (S,P):=—^- [ g(S,p)d lNls (P). 

Jt nj ^K^,J) SeC(NJ) JPoly(S) 

In other words, dfiNj is defined by putting counting measure on Cnj and by putting the 
conditional measures d^^\s (given by (9) with ip a = m^) on the 'fibers' of n. 
We then put the product measures 

dfiN,f u ...j k = dfi NJl x • • • x dfi NJk 

on the space J-N,{h,...j m ) of systems. 

We define the measure dfi^j A on Fnja anc ^ on the associated systems analogously. 
Similarly we define the measures dji^j and dji^j for the general toric and Kac-Hammersley 
ensembles, respectively. 

2.2. Expected zeros of fewnomial ensembles. We recall the probabilistic Poincare- 
Lelong formula (see for example, [SZ1, SZ3]): 

Proposition 2.1. Let (L, h) be a Hermitian line bundle on a compact Kahler manifold M . 
Let S be a subspace of H°(M,L) endowed with a Hermitian inner product and we let 7 be 
the induced Gaussian probability measure on S. Then the expected zero current of a random 
section s G S is given by 

E 7 (Z S ) = y^dd\ogU s (z : z)+ Cl (L : h) . 
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If Sj is a base-point- free linear system with Gaussian probability measure jj, for 1 < j < k 
(where 1 < k < m), then the expected value of the simultaneous zero current of k independent 
random sections Si G <Si, . . . , Sk G Sk is given by 

k 



J 7lr-,7t 



{Z su ..., Sk ) = j\ (^dd]ogUs J (z,z)+c l (L,h)\ 

.7 = 1 ^ ' 



which is a smooth form. 

Applying Proposition 2.1 to a fewnomial system S = Poly (5), we have 

Proposition 2.2. Let S 1: ...,S k be finite subsets of iVE n Z m . Then the expected zero 
current in (C*) m of k random fewnomials P± G Poly (Si), . . . , Pt G Sk is given by the smooth 
form 

^N\s 1 ,...,s k Z Pl ,...,p k = f\ i—-dd\ogU Sj (z,z) + — u FS \ . 

3=1 ^ ' 

Proof. We recall that the base point locus of a suspace S C Poly(iV) is the set of points 
at which p(z) = 0, Vp G S. Since a monomial z^ 1 ■ ■ ■ z^ 1 vanishes if and only if Zj = for 
some j such that ctj > 0, the base locus of Poly(S) is always contained in the coordinate 
hyperplances [J^ =1 {zj = 0}. Applying Proposition 2.1 to (C*) m C CF m , we obtain the 
result. □ 

COROLLARY 2.3. The expected zero current in (C*) m of a system of k random fewnomials 
of degree < N with fewnomial number f is given by 

1 k 

1 f\f-i^, „ N 



^NjZ Plt ... t P k 



C(N, f) 



where C(N,f) is given by (30). 

2.3. Mass asymptotics and fewnomial Szego kernels. We now give the asymptotics 
of the Szego kernels Hn\s We need joint asymptotics in N and S (leaving the fewnomial 
number / = | jS' j fixed). We begin with the dilated fixed spectra system (I). 

A special case of Theorem 4.1 in [SZ1] on the mass asymptotics for polynomials with 
spectra in dilates of a Newton polytope P is where P = {(3} is a single lattice point in pE. 
In this case 

U NplNm (z, z) = \m%\ 2 = Nfe- NbP ^[c + diV" 1 + c 2 N~ 2 + •••], (36) 

where 

m R I P I 2 

j=0 t> 

In (37), we can let ft be any point in the interior of pE. We also write b x — b]., for arbitrary 
(not necessarily integral) x G E: 



(37) 



b x {z) = ^2x j logx j -^Xjloglzjl 2 + log{l + \\z\\ 2 ) (x = l-\x\). (38) 

3=0 3=1 
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The first term is the symplectic potential for the Fubini-Study metric, i.e. the Legendre 
transform of the open orbit Kahler potential (see §4) We now give a precise estimate for the 
joint asymptotics of (36) using Stirling's formula. A similar analysis was done in dimension 
one in [SoZl] and in Lemma 6.2 of [SoZ3]. The kernel (36) is denoted V h N(a,z) in [SoZ2] 
and is analyzed for general toric varieties in Section 6 of that article. Since it is elementary 
we give a self-contained proof in the case of SU(m + l) polynomials (i.e., for the Fubini-Study 
metric). 

Lemma 2.4. There exists positive constants C m depending only on m such that for all a G 
(iV£) n 7L m , we have 

1 m 

log \m N a | 2 = —N b a/N + | log N - - log (f ) + R(a, N, m) , 

j=0 

where \R(a, N, m)\ < C m . 

Proof. Let x = a/N. Recalling (22), it suffices to show that 

lo § \ ^ N t™^ ( Ny ) \ = -jb( Nx i + ^j^gx j + j\ogN+R(a,N,m) (a = N- \a\) . 



Using Stirling's formula 

1 1 



n\ = V2^n n+1/2 e~ n+£n , where < e n < , (39) 

Yin + 1 Yin 



we obtain 

, ~(N + m)\(N 
log 



= 1 °S( N y lQ g( 27r ) + + log N - + 3) lo & a j + £ " ~ 

j=l j=0 j=0 



m 
~2 



log N - ZJ y*** + 2 ) log Xj + R ' 



where 



Thus 



R = M 1 + j/ N ) ~ y log(27r) +zn~Y1 £ »i 
3=1 

m 

< J2 + J) + y M 27r ) + -y2~ • 



□ 



Lemma 2.5. There exists positive constants C' m such that 

-Nb a/N + 
for allaE (#S)nZ m . 



777 

N b a/N + —\ogN-C' m < log|m^| 2 < -N b a/N + m\ogN + C' m , 
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Proof. We first suppose that a G In '■— (NT,) n Z m . The lower bound is an immediate 
consequence of Lemma 2.4. If In 7^ 0, then N > m + 1 and the maximum value of the 
convex function a4 - Sj=o 1°S (7/) 011 * s attained on the vertices of In- Thus 



1 v^, fOLi\ m, „ T 1, iV to, , T 1 



N -m 



and the upper bound follows from Lemma 2.4. 

Now suppose that a G d(NT) HZ" 1 . By a permutation of homogenous coordinates, we can 
assume without loss of generality that a = (ao, . . . , 0, . . . , 0) where Oij > 1 for < j < k. 
Let a' = (ai, . . . , a*;), z' = (zi, . . . , Zk)- We note that 

b a >/N{z') - b a/N (z) = log(l + ||*'|| 2 ) - log(l + II^H 2 ) . 

By the lower bound proved above for the monomial m^, on (C*) fe , we have 



log|m^)| 2 = log|m-(z')| 2 + log 



k 



(N + to)! 
{N + k)\ 



+ N\og 



1 + \\z' 
l + \\z\ 



/||2 



> -iV 6 a7JV (/) + - log N - C' k + (to - k) log N + N log 

= -N b a/N (z) + (m - ^ log iV - C£, 

which yields the desired lower bound when a is in the boundary of NT. 
On the other hand, by the upper bound for the monomial m^, we have 



v'l|2 



1 + \\z\ 



log \m a (z 



log|m^(z')| 2 + log 



(N + m)\ 



(N + k)\ 

< -Nb a , /N (z') + klogN + C' k + 

= -Nb a/N {z) +m log N + C' k + 
which gives the desired upper bound. 



+ iVlog 
m — k) {log N 



l + \\z'f 
1+ \\z\\ 2 
m 
N 



+ Nlog 



1 + \\z , 

1 + IUI 



'112 



to 



/cm 



N 



□ 



2.4. Proof of Theorems 1 and 2. These theorems are consequences of the following 

convergence result: 



\ogU Np \Ns(z,z) ->■ -p m.inJfe AJ/p (2;)} 



Lemma 2.6. Lei m,f,p be positive integers. Then 

1 

N ~ ^ ' - - i<j</ 

uniformly for z G (C*) m , 5 G C(p, /). 

Proof. Let 5 = {A 1 , . . . , A-^} C pT, and recall that 

/ 

n Np \Ns(z,z) = ^2\m N N p XJ (z)\ 2 . 
3=1 
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By Lemma 2.5, we have 

max{-Npb X i/ p } - C m < max{log \m" p xj \ 2 } < \ogU NplNS < max{log |m^| 2 } + log/ 
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< max{-Np b X j/N} + m log(Np) + C' m + log / 



Dividing by N, the conclusion follows. 



□ 



Remark: Lemma 2.6 is a special case of the generalization (with a stronger uniformity 
result) of Proposition 4.2 in [SZ2] to nonconvex polytopes [unpublished]. In the case where 
S is one point, an analysis of the full (i.e. not just logarithmic) asymptotics of VhN(a,z) is 
given in Section 6 of [SoZ2]. 

Proof of Theorems l-2:Let S = {A 1 , . . . , A / } C pE. By Proposition 2.2 with k = m, 



N m 'EN\NsZp u ..., Pri 



2tt 
MA 



-dd 



N 



logU Np \ NS (z,z) 



+ -o; F s 

7T 



A' 



logU Np \ NS (z,z) 



+ plog(l + \\z\\ 2 ) 



By (38) and Lemma 2.6, 



— log n NplNS (z,z) 



+ plog(l + II^H 2 ) ->■ p max (p, A p ) - (A p , log A p ) + phgp (40) 

AGS' 



uniformly, where p = (log \zi\ 2 , . . . , log |-2 m | 2 ). Theorem 2 then follows from Proposition 2.2 
and the Bedford- Taylor theorem [BT, Kl] on the continuity of the operator (ui, . . . ,Uk) h- > 
dd c u\ A • • • A dd c Uk under uniform limits. Theorem 1 is a special case of Theorem 2. □ 

Corollary3 follows immediately from Theorem 2 by averaging over the spectra in A. 

3. Zeros of random fewnomial systems: Proof of Theorem 4 

By Corollary 2.3 and the Bedford- Taylor continuity theorem for dd c u\ A • • • A dd c Uk under 
uniform limits, to prove Theorem 4, it suffices to show that 

CWT) ^ (^ogU NlSuppa {z,z)+log{l + \\z\\ 2 ) 

-> / max \(p,X J ) ~ (V,log A?)l dX 1 ■ ■ ■ dX f (41) 
y s/ j=i,...,f L J 

uniformly on compact subsets of (C*) m . 

We begin by writing the above sum as an integral. For a G A^E, we write [_a\ = 
( , • • • , [a m \ )GiVSnZ m . For a = (a 1 , ... , a / ) G (Z m n WE)', we consider the m/-cube 
of width jf 



Then 



R N , a := {(A 1 , ...,X f )e (R m ) f ■ [NX j \ = a\ 1 < j < /} . 

/ 



( ^/^ £ logn iV |s uppCT (z,^)= / log^lmf^.j^l 2 ^ 1 --.^, (42) 
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where d\ J = ml d\{ ■ ■ ■ d\ J m , and 

U N = \J {R N , a :a = {a\...,a f )e{Z m n NZ) f , a 3 ^ a 3 ' for j ^ /} . 
It then follows from (30) and the estimate Vo\(^ f AU N ) = 0(1 /N) that 



— ^— ^ogU NlSuppa (z,z) 



f 



s/ 



log 



m 



N 



(z)\ 2 dX 1 ■ ■ ■ d\ f + E N (z) , (43) 



where 



C 

\E N (z)\ < — y max 



3=1 

As in the proof of Lemma 2.6, we conclude from Lemma 2.5 that 

/ 

max{-Nbp j/N (z)} - C' m < logj^ \m%{z)\ 2 < max{- Nb p /N (z)} + m\ogN + C' m + log/. 



Therefore, there are positive constants C, C depending only on m, f such that 

\E N {z)\ < Csupb x (z) + C . 
Aes 

Lemma 3.1. Let ^ : S / x (C*) m R 6e #wen fry 

/ 

*(M = log^|mf WAJj (2 ; )| 2 . 



(44) 



(45) 



j'=i 



T/ien for all compact sets K C (C*)"\ 



1 , 
— ^(A, z) — >■ max{— 6 AJ uniformly on T,' x K . 
N j 



Proof. Let £ > be arbitrary. By (44), we can choose N such that 
1 J 

— \og^2\m^(z)\ 2 - m^{-b^ /N (z)} 



j'=i 



< £ V/3 G S / , Vz G (C*) m , ViV > iVo . 



We can choose iV large enough so we also have 

1 



| a - A| < 



iV n 



b a (z) -b x (z)\ <e Va,AGS, VzeK. 



Thus, for all (A, z) G x K and N > N , we have 



±-*(\,z) -max{-bv(z)} 
iv j 



< 



^(A, z) - max{-fe LJVA , J/Jv i 



+ 



max{-b lNX 3}/ N (z) ~ max{-6 A3 (z)} 

J 3 



< 2e. 



□ 
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The desired uniform convergence (41) follows from (38), (43), (45), and Lemma 3.1, which 
completes the proof of Theorem 4. 

The same argument gives the proof of Theoem 6. 

3.1. Computing the explicit formula: Proof of Corollary 5. For r e M m , we write 
e r = (e ri ,...,e rm ), so that = \\e p/2 \\ 2 = \e p \. 

Recalling (37), we write 

b(X;p) := b {x} (e p ) = (A, log A) - (p, A) + log (1 + | e '|) > 0. 

Therefore, 

/ max [(p,A J ) - (A J ',logA J ')l d\ l ---d\ f = log(l + |e"|) - / min 6( A J ; ; p) dX 1 ■ ■ ■ d\ f . 
J-zf j=h-J L J Jxf j=i,...,f 

(46) 

We shall use the following elementary probability formula: Let X be a non-negative ran- 
dom variable on a probability space (Jl,dP), and let D x (t) : = P(X < t) be its distribution 
function. The expected value of X is given by 

//*00 PT 
XdP= / tdD x (t)=]im / tdD x (t), 
Jo r ^°° Jo 

where 

/ tdD x (t)=rD x (r)- [ D x (t)dt= [ [D x (r) - D x (t)} dt . 
Jo Jo Jo 

Letting r — > oo, we have by Lebesgue monotone convergence 



E 

We let 



XdP = J [1-D x (t)]dt. (47) 



D b (t;p) :=P{AGS:ft(A;p) < t} 



be the distribution function for b(-;p), where dP(X) = m\dX 1 ---dX m . The distribution 
function for the random variable 

X(X 1 , ...,X f ) -min^A 1 ),...,^)} 

on (with the product measure dP^X 1 ) ■ ■ -dP(X^)) is given by 

£>* = 1 - (1 - £> 6 )'. 

It then follows from (46)-(47) that 

/ max \{p,X J ) ~ (X j ,\ogX j )] dX 1 ■ ■ ■ dX f = log(l + \e p \) - [ [1 - D b (t; p)] f dt . (48) 
Corollary 5 follows immediately from Theorem 4 and (48). □ 
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3.1.1. The dimension 1 case. We now further evaluate A, when the dimension m — 1. In 
this case, 

b(X;p) = A log A + (1 - A) log(l — X) — p\ + log(l + e p ) , < A < 1, peR. 

Since b is a convex function of A (taking the minimum value when A = e p /(l + e p ) ), we 
have Db(t; p) = g(t, p) — g(t, p) for t > 0, where g(-, p) < g(-, p) are the branches of b(-, p) _1 . 
Precisely, g = g(t, p), g = g(t, p) are given by 

< g < g< 1, 

b(g;p)=t if t < log(l + e p ) , b(g; p) = if t > log(l + e p ) , 
6(5^; p) = t if t < log(l + e-P) , 6(5^; p) = 1 if t > log(l + e"") . 

We have the symmetry b(X; p) = 6(1 — A; — p), and hence g(t, p) = 1 — (?(t, — p). Therefore, 

D b (t;p) = l-g(t,p)-g{t,-p), (49) 

where g(-, p) : [0, +00) — > [0, e p / (1 + e p )] is given by: 

b(g(t,p),p)=t, if 0<t<log(l + e p ), 

g(t,p) = 0, if t>log(l + e p ). ^ 

4. General toric Kahler potentials 

We now sketch the proof of Theorem 7. It is almost the same as in the Fubini-Study case 
but requires the generalization of Lemma 2.5 and then Lemma 2.6. 

As discussed in [SoZ2], the toric norming constants can be written in terms of the sym- 
plectic potential as follows: 

Qg„( V 4V v )((x) = jf e -"(M*)+<#-*,iW(*)><fe. (51) 

Here, p ip {e p ^ 2 ) = V p <p(e P//2 ) is the moment map determined by if. Applying steepest descent 
to the integral, we find that there exists only one critical point at x — p ip (e p ^ 2 ), and we 
conclude that 

logQ Gjv( ^)(a) = iv (^) + O (^jP) (52) 

uniformly [SoZ3, (25)]. 

The logarithmic asymptotics (52) is the only non-obvious aspect of the logarithmic mass 
asymptotics. The Szego kernel for a single lattice point (on the diagonal) equals 



e (a,p) e -AV(e p/2 ) 

f-N,Q\ a Ke"'-,e"'-) = — . 

WG N (<p,dV v ) 



U NiQla (e p / 2 ,e p / 2 ) 



The analogue of Lemma 2.5 for a general Kahler potential is 

\ogU N , QlNx (e p / 2 ,e p / 2 ) = N((x,p)- V (e p ' 2 ) - u v (x)) + O(logiV) , (53) 

which follows from (52) and [SoZ3, (55)]. 

For a fewnomial Szego kernel with a finite set S of lattice points, the analogue of Lemma 
2.6 is that 

^ log U NiQlNS (e p f 2 , e p l 2 ) = max((A, p) - u^X)) - V {e p l 2 ) + O (^jpj . (54) 
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The proof is the same as that of Lemma 2.6, using (53). With this modification, the remain- 
der of the proof of Theorem 7 is the same as that of Theorem 4. 

5. The T m ensemble 

Finally, we indicate the modifications needed to deal with the fewnomial Kac-Hammersley 
ensemble (10) . This is quite different from the case of pluri-subharmonic weights because 
the Szego kernel has quite different (much weaker) asymptotic properties. But for fewnomial 
Szego kernels the distinction is not too severe. 

In this case, we use the L 2 norm || • || T m on the real torus rather than the Fubini-Study 
norm. We therefore have 

E Tm{s (\p( z )\ 2 Tm ) = J2 E(Kh)xMxi&) ■ 

Since E(A Q ,A (3 ) = 5^, we have: 

E Tm]s (\P(z)\ 2 Tm ) =J2 \Xa(z)\ 2 = U Tmls (z,z) , (55) 

where n T m| S is the orthogonal projection onto Poly(S') C L 2 ((C*) m , 5 T m). It then follows 
by expressing the Gaussian in spherical coordinates that the expectation in the fewnomial 
Kac-Hammersley ensemble is given by 

E KH (\P(z)\ 2 Tm ) = ^E KH (\P(z)\ 2 Tm ) = ^U Tmls (z,z) . 

It is clear that 

U Tmls (z,w) = J2( z ^) a - (56) 

aeS 

Therefore, 

The potential in this case is 

F ^ Z):= Cihn ^ ^gU T m, s (z,z). (58) 



Proposition 5.1. 



lim -J- F f N (e p/2 ) = [ max{(ii, p) . . . , (x f , p)}dx 1 --- dx f . 



Outline of the proof: Indeed, 

log Y,e {p ' a) = N\ogJ2 e{p,a/N) ~ Nm&x{(p,a/N)}. (59) 



Hence, 

F N^ P ) = CiWT) E log^e^-Tv/" m a x{(x l ,p)...,(xf,p)}dx l ---dxf. 



□ 
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We note that for each (x 1 , . . . ,x^), the function M( x i ) ... jX m)(^) = max-^x 1 , /o) . . . , (x? , p)} 
is a piecewise linear convex function. It follows that the integral defines a convex function 
of p. 

In dimension one, if all x^ > 0, 



max-fpx 1 , . . . , px* 

Hence, 



/jmaxji 1 ,...,^}, p > 0, 
pmin-fx 1 , . . . , x?}, p < 0. 



Np{ Jj ^ maxja; 1 , . . . , x^dx 1 ■ ■ ■ dx?}, p > 

J$(e") ~ { ' (60) 

iVp{ Jj ± y minja; 1 , . . . , x-^jjdx 1 ■ ■ ■ dx?}, p < 

Thus, F^(e p ) is piecewise linear in p with a corner at p = 0. In dimension one, 



i E^Zp*) = ^ddF* 5 51 . (61) 
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